Junction between surfaces of two topological insulators 
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We study the properties of a line junction which separates the surfaces of two three-dimensional 
topological insulators. The velocities of the Dirac electrons on the two surfaces may be unequal and 
may even have opposite signs. For a time reversal invariant system, we show that the line junction is 
characterized by an arbitrary parameter a which determines the scattering from the junction. If the 
surface velocities have the same sign, we show that there can be edge states which propagate along 
the line junction with a velocity and orientation of the spin which depend on a and the ratio of the 
velocities. Next, we study what happens if the two surfaces are at an angle (j> with respect to each 
other. We study the scattering and differential conductance through the line junction as functions 
of (j) and a. We also find that there are edge states which propagate along the line junction with 
a velocity and spin orientation which depend on 0. Finally, if the surface velocities have opposite 
signs, we find that the electrons must transmit into the two-dimensional interface separating the 
two topological insulators. 

PACS numbers: 73.20.-r,73.40.-c 
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I. INTRODUCTION 

Recent years have witnessed extensive theoretical^^^ 
and experimental^"— studies of a class of two-and three- 
dimensional materials called topological insulators; for 
reviews, see Refs. [H and [IS- A topological insula- 
tor (TI) is a material which is gapped in the bulk but 
has gapless states at the surface which is one- or two- 
dimensional if the TI is two- or three-dimensional. Fur- 
ther, the electrons have strong spin-orbit coupling, and 
the surface states are described by a massless Dirac equa- 
tion (in one or two dimensions) in which the directions 
of the spin angular momentum and linear momentum 
are tied to each other; they both lie in the plane of the 
surface and they are perpendicular to each other. The 
three-dimensional TIs come in two classes, strong and 
weak, which respectively have an odd and even number 
of massless Dirac cones at the surface^"— i^. In strong 
TIs, the gaplessness of an odd number of Dirac cones is 
protected by time reversal invariance; time reversal in- 
variant perturbations such as non-magnetic disorder and 
lattice imperfections do not produce a gap, while time re- 
versal breaking terms such as magnetic impurities or an 
external magnetic field can produce a gap for the surface 
states. For materials such as HgTe, Bi2Se3 and Bi2Te3, 
specific surfaces have been found which have a single 
Dirac cone at one point of the two-dimensional Brillouin 
zone^ d^d^d'^' J^. Several features of these surface Dirac 
electrons have been investigated recently. These include 
the existence of Majorana fermion modes at a magnet- 
superconductor interface22ri2^, anomalous magnetoresis- 
tance of ferromagnet-ferromagnet junctions^!, realiza- 
tion of a switch by magnetically tuning the transport of 
the electrons by a proximate ferromagnetic film^'^, spin 
textures with chiral properties'^"—, and the anisotropy 
of spin polarized scanning tunneling microscope (STM) 
tunneling into the surface of a TI^^. Some general proper- 
ties of the surface states and their scattering at the edges 
of two- and three-dimensional TIs have been studied in 



Ref. 121 

A recent paper^ has studied what happens when two 
three-dimensional TIs are brought close to each other so 
that their surfaces touch each other along a line; we will 
call this a line junction. This produces a system in which 
two two-dimensional surfaces governed by massless Dirac 
equations, with possibly different velocities vi and V2, 
share a one-dimensional boundary; similar studies have 
been carried out in Refs. |29| andlSOl and in one dimension 
in Ref. l3ll . If a plane wave is incident on the boundary 
from one side with an angle of incidence ^, it will be 
reflected and transmitted with certain amplitudes which 
depend on 9 and the ratio of velocities vxjvi. Ref. [23 
made the interesting observation that if the velocities v\ 
and f 2 have the same sign, the scattering problem can be 
solved in a straightforward way, but if they have opposite 
velocities, there is a difficulty; at least for the case of 
normal incidence, conservation of spin implies that the 
refiection and transmission amplitudes must both vanish. 
The way out of this difficulty is to consider two more 
sets of surface states which lie at the interface of the 
two TIs; then these states necessarily have some gapless 
states, even if one takes into account tunneling between 
them. Ref. |28| proposed that a low-energy plane wave 
which is incident normally on the line junction must get 
transmitted into these interface states. 

In this paper, we will generalize the study of a junction 
between two TI surfaces in several ways. In Sec. II, we 
will consider the case of velocities with the same sign. We 
find that the most general condition at the line junction 
which conserves the current has one arbitrary parameter 
a if we assume time reversal invariance. (The condition 
assumed in Ref. [28 corresponds to the special case of 
a = 0). We will compute the refiection and transmission 
amplitudes across the line junction as functions of ^, a 
and the ratio of the velocities on the two sides of the line 
junction. We will also show that depending on the values 
of a and the ratio of the two velocities, there can be edge 
states which travel along the line junction with a momen- 



turn k but decay exponentially perpendicular to it. The 
velocity of the edge states along the junction is smaller 
than the velocities of the plane waves propagating far 
from the junction, while the spin orientation of the edge 
states has a component pointing out of the surface in con- 
trast to the plane waves whose spins lie in the surface. 
In Sec. Ill, we will consider a TI which has two surfaces 
which meet at a line junction at an angle 0; we will show 
that this provides an interesting scattering problem in 
which the reflection and transmission amplitudes depend 
on 0, the bending angle and the parameter a. We will 
use the transmission probability to compute the differen- 
tial conductance across the line junction. Further, we will 
show that for 7^ 0, there are edge states which travel 
along the line junction with a momentum k but decay 
exponentially away from the junction; the velocity and 
the spin orientation of these states depend on (j). In Sec. 
IV, we will examine the case of two TI surfaces where the 
velocities have opposite signs, and we will introduce an 
interface with two additional surfaces. We will explicitly 
derive the energy spectrum and wave functions at the in- 
terface assuming a simple form of the tunneling between 
the interface surfaces. We will then demonstrate that a 
wave incident on the surface indeed transmits into these 
states although the reflection and the transmission am- 
plitudes on the surface are generally non-zero if the angle 
of incidence 9^0. In Sec. V, we will point out some 
ways by which the edge states moving along the junction 
may be experimentally detected, and we will make some 
concluding remarks. In the Appendix, we will study the 
most general boundary condition which is consistent with 
the Hermiticity of the Hamiltonian and the conservation 
of the current at the junction in the model discussed in 
Sec. II. We will show that the general condition involves 
four real parameters. However, if we assume that the 
system is invariant under time reversal, then there can 
be only one parameter; physically this corresponds to the 
strength of a potential barrier which may be present at 
the junction. 




FIG. 1: (Color online) Schematic picture of two regions of the 
X — y plane, labeled 1 and 2, which are separated by a line 
junction at y = where a thin barrier may be present. A wave 
incident (I) from region 1 and the corresponding reflected {R) 
and transmitted (T) waves are shown. 



Eq. ([IJ act on wave functions -04 1 where i = 1 and 2 for 
2/ < and > 0, respectively. (We will generally set ft = 1 
everywhere) . 

Hamiltonians defined in disjoint regions, such as the 
ones given in Eq. ([T]), do not define a system completely; 
this has been emphasized earlier for the Dirac equation in 
one dimension"^ ^. To define the system fully, we need to 
specify the boundary conditions that the wave functions 
must satisfy at y = in order to ensure that the Hamil- 
tonian H of the entire system, given by Hi for y < and 
H2 for J/ > 0, is Hermitian, namely, that 
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{Hx')^X (2) 



for any two wave functions x ^^id x' which satisfy the 
boundary conditions. On doing an integration by parts 
with dy, we find that Eq. ^ requires 



VI ((x'i)V-xi) 



y^o- 



V2 ((xi)t^"X2) 



(3) 



II. JUNCTION BETWEEN SURFACE 
VELOCITIES WITH SAME SIGN 

In this section, we will study a system with a line 
junction which separates the surfaces of two three- 
dimensional TIs, labeled 1 and 2. The surfaces define 
the X — y plane, and the line junction lies along y = 0, as 
shown in Fig.[TJ a thin barrier may be present along y = 
as we will discuss below. The Hamiltonians in the two re- 
gions, y < and y > 0, are given by the two-dimensional 
massless Dirac form 



Hi = -ivi {a^dy ~ CT^Sa;), 
H2 = -iV2 [a'^dy - a^d,), 



(1) 



where a"" denote Pauli matrices, and vi, V2 denote the 
Fermi velocities on the two surfaces respectively. We will 
assume here that ui,W2 > 0. The Hamiltonians Hi in 



We may now ask: what is the most general linear 
relation between (xi)y^.o- ((Xi)y^o-) and (x2)y^o+ 
{{X2)y^o+) which satisfies Eq. ([3])? The answer to this is 
discussed in detail in the Appendix; we show there that 
the general linear relation involves four arbitrary real pa- 
rameters, and we provide a physical understanding of 
these parameters. We then argue that if the system is 
invariant under time reversal, the linear relation involves 
only one real parameter a and is given by 



(X2) 



y^o+ 



(Xi) 



a->o- 



(4) 



and similarly between (x2)a-^o+ and (xi)a-+o-- We note 
that changing a ^ a + n has no effect on any physical 
quantities (such as reflection and transmission probabili- 
ties) since this is just equivalent to changing 7/12 ~^ ~V'2- 
In the rest of this paper, we will assume that a lies in the 



range -7r/2 < a < 7r/2. In Refs. l28|-|3l|, the relation be- 
tween the wave functions at the hne junction was taken 
to be 



iX2)y- 



>0+ 



(Xi) 



y^Q-- 



(5) 



Clearly this is a special case of Eq. (|3]) corresponding to 
a = 0. 

Given the form of the Hamiltonians in Eq. ([T]), the 
current density Ji = {Jix^Jiy) on surface i {— 1,2) is 
given by 



J^ 



V, {~4a^i,,,i,\a^^i). 



(6) 



[This can be derived by using the equation of motion 
'idt'^i = Hiii^i and the equation of continuity dtPi + V • 
Ji = 0, where the charge density is pi = i'ji'i-] Current 
conservation in the y direction at the line junction at 
y — implies that 



Vi {■lpla^^pi)y^o- = V2 (V'|cr^V'2)iy->0+, 



(7) 



We now see that this is ensured by the same condition 
^ which ensures Hermiticity of the Hamiltonian H. 

The system described by Eqs. (P) and dU is invariant 
under two discrete symmetries, namely, parity and time 
reversal. Under parity, x — >■ —x and ip — >■ a^tp- Under 
time reversal, t —)' —t and ip —> a^tp* . 



A. Scattering 

For an electron far from the line junction, the eigen- 
states of the Hamiltonian are plane waves labeled by a 
momentum k ~ {kx,ky). For a velocity Vi > 0, the en- 
ergy spectrum and wave functions are given by 






J{k^x+kyy-E+t) 



^i(k^x+kyy — E-i) 



(8) 



for E^ ^ Vik and E^ = —Vik respectively, where cos 9 — 



ky/k, sin 6' = k^/k, and k — Jk"^ + ky. Note that the 

directions of the spin and momentum are tied to each 
other; for positive energy states, they are related by ct = 
fc X z, where x, y and z form a right-handed coordinate 
frame. 

Now let a plane wave be incident on the line junction 
from the left {y < 0) with momentum fc/ = (kix, kjy) and 
positive energy _E_|_ — vi\ki\. The reflected and trans- 
mitted waves will have momenta kn = {knx^kny) and 
kx = (kTxjkry) and amplitudes r and t respectively. We 
can determine all these quantities in terms of the incident 
momentum kj using conservation of energy, momentum 
in the x direction (due to translation invariance in that 



direction) , and current conservation in the y direction at 
the line junction. Conservation of energy implies that 

Vl^kJ^+kJy = Vl^kl^+kly = V2^kl^+kly, whUC 

conservation of momentum in the x direction implies that 

kix = kRx = kTx- (9) 

These equations imply that 

k-Ry ^ rC/y, 



kry = — 



L^vl{kj^ + kjy) -vikl 



(10) 



assuming that the quantity inside the square root in the 
second equation is positive; otherwise k^y will be imagi- 
nary and we will have total internal reflection on the left 
side of the line junction^^. Thus kj^ and fc^ are fixed 
in terms of kj. If we define 9 = ta.n^^ {kix/kjy) and 
e' = ta.n-^{kTx/kTy), where -7r/2 < 9, 9' < 7r/2, we 
have 



1 



Vl 



smf 



1 



V2 



(11) 



Finally, current conservation in the y direction at the line 
junction at y = implies that 



vikiy (1 - |rp) = V2kTy \t\'^ 



(12) 



The reflection and transmission amplitudes r and t 
can be obtained for the general boundary condition in 
Eq. ((H). The general expressions for r and t are 



t= J^ 



cos 9 



9'-e)/2 



V2 cosa cos(^^) -I- isina cos(^-2^)' 
g sina sin(^i^) + icosa sin(^-^) 
cosa cos(^i^) + isina cos(^^) 



(13) 



For equal velocities vi = V2, we have 9' — 9 and Eq. P^ 
simplifies to 



t = 



cos 6* 



cos a cos a + I sm a 
e~'^ sin a sin 



cos a cos ( 



I sm a 



(14) 



Given the transmission amplitude i, we can compute 
the conductance across the line junction. This will be 
discussed in Sec. IV for a more general model in which 
the surfaces 1 and 2 are at an angle with respect to each 
other. 



B. Edge states 

It turns out that for certain ranges of values of Wi/u2 
and a, there are edge states which propagate along the 
line junction (i.e., along the x direction) with a momen- 
tum k and whose wave functions decay exponentially as 



y — >■ ±00. At the line junction, let us look for unnormal- 
ized wave functions of the form 



^1 



1 

7i 

72 
73 



J{kx-Et)+xiy 



in region 1, 



^2 = [ J] Skx^Et)'X2V in j.ggion 2, (15) 



where 71, 72, 73 can be complex numbers but the Xi 
must be re al and p ositive. The energy is given by 
E = ±vi^k'^ — Xi — ±V2\/k^ — xh ^^ must have 
Xi) X2 ^ 1^1- We now demand that the wave functions in 
Eq. (|15l) be eigenstates of Eq. ([T]) and that they satisfy 
Eq. Q; this gives us the values of the 7^ and Xi- Let 
us parametrize Xi/^ = sin^i and X2/fc = sin02, where 
— 7r/2 < 9i, 62 < tt/2. We further assume that vi < V2, 
and define (3 — cos~^(wi/w2), where < /3 < 7r/2. 
We then find that edge states exist under the follow- 
ing conditions. For f3/2 < a < tt/2, the energy is 
negative and given by -E = —vw/k"^ — xf , while for 
— 7r/2 < a < — /?/2, the energy is positive and given 
by i!^ = vi^k'^ — xl- In both cases, the magnitude of 
the velocity, \E/k\, is given by 



V = 



wi?;2| sin(2Q!)| 



vw 



— 2wiW2 cos(2q;) 



(16) 



Further, 



tan 01 
tan02 



sgn{E/k) 



vi — V2 cos(2a) 
V2 sin(2a) 



fjp/i,\ ^^2 - wi cos(2a) 

sgn{E k) —-— — . (17 

Vi sm(2a) 



where sgn denotes the signum function: sgn{z) = +1 if 
z > and — 1 if z < 0. Note that there are no edge states 
if -13/2 <a< (3/2. 

We can use the above expressions to show that the 
velocity v is smaller than both vi and V2. To state this 
differently, for a given value of the a:-momcntum fc, the 
energy of the edge states has a smaller magnitude than 
the energy of the plane waves in both regions 1 and 2. 
(At the two points a — ±/3/2, we find that v — vi). 

The various expressions above simplify if vi = V2- 
Then there are edge states in the entire range — 7r/2 < 
a < 7r/2, with v = V2\ sina| and 9i = 02 = —sgn{E/k) a. 
In Fig. m we present a plot of the velocity of the edge 
states as a function of a, in the range — 7r/2 < a < 7r/2, 
for Ui/w2 — 1 and 0.5. 

We emphasize that the sign of the energy E of the edge 
states depends only on the sign of a, and not on the sign 
of the momentum k. This is consistent with the time 
reversal invariance of the system; under time reversal, 
the momentum reverses sign but the energy remains the 
same. 

After obtaining the values of the 7j in Eq. (IT5|) . we 
can find the orientation of the spin of the edge states. 
For a spin- 1/2 particle whose spin points along a unit 




FIG. 2: (Color online) Plot of v/vi versus a for v\/v2 = 1 
(solid) and 0.5 (dashed). 



vector h = (sin a cos 6, sin a sin 6, cos a) , we know that the 
normalized wave function is given by the transpose of 
(cos(a/2),sin(a/2)e*''). Using this we find that the spin 
of the edge states points along the direction — sin 6'i z =F 
cos Q\ y in region 1 and along the direction sin O2 z ^ 
cos 02 y in region 2, where the =F signs refer to ii^ > and 
< respectively. Thus the spin of the edge states lies in 
the y — z plane and not in the x — y plane which is the 
surface of the TIs. 

To conclude this section, we see that several properties 
of the edge states depend only on the dimensionless pa- 
rameters a and v\/v2- This is true of the velocity E/k, 
the spin orientations in the two regions, and the product 
of k and the exponential decay lengths 1/xi Etnd l/x2. 



III. JUNCTION BETWEEN SURFACES AT AN 
ANGLE 

We will now study what happens when there is a line 
junction between two surfaces which are boundaries of 
the same TI, but the surfaces are at an angle with re- 
spect to each other; this is shown in Fig. 121 This is a real- 
istic situation since any finite-sized TI must be bounded 
by surfaces meeting along lines. We will assume that the 
velocity v is the same and is positive on both surfaces. 
Region 1 forms the x — y plane, while region 2 forms the 
x — y' plane where y' = cos 4> y — sin (jj z. In regions 1 and 
2, y < and y' > respectively; the line junction lies at 
y = y' — 0, and we will again allow for a thin barrier to 
be present there. The Hamiltonians in the two regions 
are given by 



Hi = -iv [a'^dy - a^dx), 
H2 = -iv [cr^'dy' - a^'dx) 



(18) 






y 
y' 



X 




FIG. 3: (Color online) Schematic picture of two regions, la- 
beled 1 and 2, separated by a line junction a.t y = y = Q 
where a thin barrier may be present; the unit vectors nor- 
mal to the two regions are at an angle 4> with respect to each 
other. A wave incident (/) from region 1 and reflected [R) 
and transmitted (T) waves are shown. 



where a'^ = cos (/) ct^ — sin (j) a^ . This model is of physical 
interest since any three-dimensional TI of finite size must 
necessarily have surfaces meeting along lines; our model 
describes the region around one such line. We will take 
<f> to lie in the range — tt < (j) < tt, although real systems 
may only allow a more restricted range of (f). The system 
described in Eq. (I18|) is invariant under the parity and 
time reversal transformations. 

We may now ask what the general linear condition 
must be between {ipi)y^a- and {ip2)y'^o+ so that the 
Hamiltonian of the complete system is Hermitian. Fol- 
lowing an analysis similar to the one given in the Ap- 
pendix for the model in Sec. II, we find that the gen- 
eral relation involves four parameters. But if we de- 
mand that the system be invariant under time rever- 
sal, we find that only one parameter can appear, namely, 
{tp2)y'-^o+ = e^^°"^ (V'i)a^o-- This also ensures current 
conservation at the line junction, namely. 



(V'Ia^Vl)y^O- = (V'lff''V'2) 



a'->o+- 



(19) 



A. Scattering 

For plane waves, the energy spectrum and wave func- 
tions in region 1 have the form given in Eq. (|S]), while in 
region 2, we obtain 

, _ e"'*^ /^ cos 6* -H I sin 6* cos \ i(k^x+kyv'-E+t) 



V'- 



V2 



1 — sin 9 sin ( 

cos 6* + i sin cos i 
— 1 — sin 6* sin (i) 



^i(k^x+kyy' —E-t) 



(20) 



where E± = ±v^/k'^ + k'l and 9 = tan ^{kx/ky). With 



the normalization given in Eq. (|20p . the currents in the 
y' direction ; 
and V'lcr^i/'- 



y' direction are given by ip^a^ip^ — cos 9 (1 — sin 9 sin ( 



We can then solve a scattering problem in which a wave 
is incident on the line junction from region 1 with momen- 
tum ki = {kix,kiy), and gets reflected and transmitted 
with amplitudes r and t respectively. These amplitudes 
depend on the incident angle 0, the bending angle cf), and 
the parameter a. The general expression for t is 

2 cos 6* 





Ai cos a + iA2 sin a ' 


Ai -- 


= cos 9 + isiii9cos(j) + e^**(l — sin 6* sin 0), 


A2 -- 


= 1 -sin6'sin0 + e"*''(cos6' + isin6'cos(/)).(21) 



We note that the transmission probability |ip remains 
the same if we simultaneously change 9 -^ —9, a -^ —a 
and (p —>■ —(p. As a result, the conductance (to be dis- 
cussed in the next section) remains the same if we simul- 
taneously change a -^ —a and (f) ^>- —cf), but generally 
changes if a — > — a or — > — separately. 



B. Conductance 

We will now compute the conductance across the line 
junction. Suppose that the system is coupled at the far 
left {y — > — cx)) and far right {y' — ?► cx)) to reservoirs 
with chemical potentials fiL and fij^ respectively; we will 
assume that the reservoirs are at zero temperature. Elec- 
trons coming from one reservoir transmit across the line 
junction and go to the other reservoir with a probability 
given by |ip given in Eq. (HI]). Assuming that the system 
has a large width in the x direction given by W, the net 
current going from left to right is given by 



/ = eW 



dk rik 

=;^ |tp t;cos0 (1 - sine sin 0), (22) 
(27r)^ 



where e is the charge of an electron, and the factor of 
V cos 0(1 — sin 9 sin 0) appears because we are interested in 
the y' component of the transmitted current in region 2. 

For the integral in Eq. (|22|) . the energy E ~ hv^/k'^ + ky 

goes from fiji to fij^, assuming that fiji < fij^, and the 
angle 9 = tan"^ (kx/ky) goes from — 7r/2 to 7r/2. The 
voltage applied to a reservoir is related to its chemical 
potential as /i = qV; in the zero bias limit, fiL, M-R ~^ M: 
we obtain the differential conductance 



G 



dl 
dV 



v(27rft)2 



7r/2 



t/2 



d9 |t|^ 



(1- 



(23) 

Using Eqs. ([21]) and ([23| . we can calculate G as a function 
of a and (p. For a — (f> — 0, we have t = 1 and we obtain 
the conductance Go — 2e^W fi/ {v{2'Kh)^) . 

Fig. 0] shows a plot of the conductance (in units of Go) 
as a function of a, in the range —tt/2 < a < 7r/2, for 
(f) — (sohd) and tt/2 (dashed). Fig. [S] shows a plot of 
the conductance as a function of (p, in the range —tt/2 < 
<P ^ ^/2, for a = 0, 7r/4 and 7r/2. Interestingly, we find 
that the conductance is unity (in units of Go) on the line 
a = — </'/2, where — 7r/2 < a < tt/2. 




FIG. 4: (Color online) Plot of G/Gq versus a for (/> = (solid) 
and 7r/2 (dashed). 




where 7 can be complex but x must be real and positive; 
the above wave functions satisfy the boundary condition 
given above. Demanding that these wave functions be 
eigenstates of Eq. ([T5| , we find that there is an edge state 
for any non-zero values of k and <j). Since — tt < </> < tt 
and 7^ 0, we can use the facts that cos((/)/2) is always 
positive and that sin(0/2) and (f> have the same sign. We 
then find that if > 0, the edge state energy is negative 
and given by E 



'V\/W- — x^7 while if < 0, the energy 
is positive and given by E 
we obtain the relation 



v^Jk?' — x^ ■ In either case, 



yjk"^ — x^ -^/fc^ — x^ + ^1 sin( 



k-x 
whose solution is given by 



k cos <^ + X 



X = |fc sin((/./2)|. 



(25) 



(26) 



For instance, x = |fc'/'/2| for (/> ^- and = |fc|/-\/2 for 
(\) — ±7r/2. The quantity 7 in Eq. (1241) is given by 



7 



cos(0/2) 



sgn(k(f)) — sin(0/2) 



(27) 



We conclude that several properties of the edge states 
depend only on the bending angle (p and on the sign of 
the momentum k. The product of the exponential de- 
cay length 1/x and |A;| is equal to 1/| sin(0/2)|. The 
magnitude of the velocity, \E/k\, is given by vcos{(l)/2). 
We may also consider the direction of the spin for the 
wave functions in Eq. (|24|) with 7 given by Eq. (j27|) . 
We find that the spin points along the direction h — 
cos{(j)/2) y — sm{(j)/2) z if fc(/) > and along ~n — 
— cos{(j)/2) y + sm{(j)/2) z if k(j> < 0. It is interest- 
ing to note that n lies half-way between y and y' — 
cos <j) y — svncj) z. 



IV. JUNCTION BETWEEN SURFACE 
VELOCITIES WITH OPPOSITE SIGNS 



FIG. 5: (Color online) Plot of G/Go versus for a = (solid), 
7r/4 (dashed) and 7r/2 (dotted). 



C. Edge states 

For any non-zero value of 0, we find that there are edge 
states which propagate along the line junction with a mo- 
mentum k and whose wave functions decay exponentially 
as either y — > — cx) or y' — >■ 00. Let us set a = 0, so that 
we have the boundary condition {'ip2)y'^Q+ — (V'i)j,^o- 
at the line junction. We now look for unnormalized wave 
functions of the form 



xpi 



J{kx-Et)+xy 



^2 = ( M e'^'^^'-^'^-^'y' 



in region 1, 



in region 2, (24) 



Let us now study the case where there are two regions 
as in Sec. |TT1 but the velocities vi and V2 in those two 
regions have opposite signs^S. To be specific, let us as- 
sume that 1^1 > and V2 < 0, and we send in a wave from 
region 1 with positive energy as before. Let us consider 
normal incidence, so that kjx = knx = kxx = 0. We then 
see that the operator cr^ commutes with both Hi and 7?2- 
For kiy > and positive energy E = vikjy, the incident 
wave function is an eigenstate of a^ with eigenvalue 1; 
hence the reflected and transmitted waves must also be 
eigenstates of a^ with eigenvalue 1. However, this im- 
plies that the reflected wave must have kf/y — kjy and the 
transmitted wave must have kxy — {vi/v2)kiy; neither of 
this is possible since a reflected wave must have kj^y < 
and a transmitted wave must have kxy > 0. We there- 
fore see that the reflected and transmitted amplitudes 
must both vanish, leading to a difficulty in conserving 
the probability. 



■y = z = 






y 



X 




FIG. 6: (Color online) Schematic picture of two regions of the 
X — y plane, labeled 1 and 2, separated by a line junction at 
y = z = 0, and two additional surfaces, labeled 3 and 4, which 
lie in the x — z plane and are at the interface of the two TIs. 
A wave incident (7) from region 1 at 61 = and transmitted 
waves in regions 3 and 4 are shown. 



One way out of this difficulty is to look for a linear 



relation between (ipi) 



y^o- 



and (-0: 



2)y-^Q+ 



which satisfies 



current conservation as in Eq. ([7|). We find that the re- 
lation 



(0: 



'2)y^Q+ 



)y^0- 



(28) 



works. However, there does not seem to be a physical 
way of justifying the factor of a^ on the right hand side 
of Eq. ([28)) . although the factor of e"*"*^"" can be justified 
in the same way as we did in the Appendix for the case 
of velocities with the same sign. In fact, the factor of a^ 
in Eq. (|28l) makes it non-invariant under the transforma- 
tion ipi^2 — >■ cr^^i_2 which was used above to show that 
the reflection and transmission amplitudes must vanish 
for normal incidence. The factor of a^ would also break 
invariance under time reversal. 

An alternative and physically appealing way out of the 
difficulty in conserving probability was proposed in Ref. 
[28, by realizing that in the case of opposite velocities, the 
waves must transmit into the interface regions 3 and 4 
separating the two TIs, as depicted in Fig. [51 Regions 3 
and 1 will both be taken to be surfaces of the first TI; we 
therefore assume that the velocity in both these regions 
is equal to vi . Similarly, regions 4 and 2 are both surfaces 
of the second TI and therefore have the same velocity V2 ■ 
The line junction is therefore a boundary of four surfaces, 
namely, 1, 2, 3 and 4; the location of the line junction 
is given hy y = for surfaces 1 and 2, and z = for 
surfaces 3 and 4 (with z becoming negative as one moves 
away from the line junction on these two surfaces). As we 
will discuss below, it will be convenient to assume that 
the line junction has a small width separating the first 
and second TIs. To see how a wave incident from region 
1 can transmit into regions 3 and 4, we must first show 
that the interfaces 3 and 4 must necessarily have gapless 
states; otherwise energy conservation would not allow a 
very low-energy incident wave in region 1 to transmit 
into the interfaces 3 and 4. Let us therefore consider the 
Hamiltonian describing an electron moving on surfaces 



3 and 4. Since 3 and 4 both lie in the x — z plane, 
but the unit vectors normal to those surfaces point in 
the y and —y directions respectively, the corresponding 
Hamiltonians are given by 



Hi = iv2 {(T'd,j, - a'-'d:,). 



(29) 



Next, we examine what boundary condition must be 
imposed at the line junction so that the complete Hamil- 
tonian is Hermitian. Since the system consists of four 
Hamiltonians, ifi, H2, H^ and H4, the general boundary 
condition which ensures Hermiticity will be much more 
involved than the one discussed in the Appendix for the 
case of two Hamiltonians. We will therefore consider the 
simplest possible boundary condition. Let us suppose 
that in the region of the line junction, the first TI (with 
surfaces 1 and 3) is separated from the second TI (with 
surfaces 2 and 4) by some distance; hence ip^ will be con- 
nected to only ipi, and -04 will be connected to only ■02- 
Let us further assume that there is no potential barrier 
between surfaces 1 and 3 and between surfaces 2 and 4, 
so that the parameter a introduced in Sec. II is now 
equal to zero. We then obtain the simple conditions 



ips = ■01 and ^04 = V'2 



(30) 



at the line junction. (Note that this gives four equa- 
tions since the ipi are all two-component spinors). This 
condition automatically implies that the incoming and 
outgoing currents normal to the line junction are equal 
to each other, namely, 

vi'ipla'^'ipi - U20|cr''02 = viipla'^ipa - V2ipl(j''ipA- (31) 

Finally, we introduce tunneling between surfaces 3 and 
4. Assuming rotational invariance in the x— z plane, the 
equations of motion take the forni2^ 

idtips == ~ivi {a'-d^ - a''dz)ip3 + gip4, 
idtipA = tV2 [o'dx - a"a^)V'4 + 5^3, (32) 

where g denotes the tunneling amplitude; we will assume 
that g is real and positive. In defining our model, we 
have made the conceptually simple assumption that the 
region of the line junction is distinct from the tunneling 
region so that the tunneling does not modiiy the bound- 
ary condition introduced in Eq. (1501) . (We note that the 
system described in Eq. (1321) is invariant under the parity 
and time reversal transformations defined in Sec. HIl) . 

To make further analytical progress, let us assume that 
V2 — ~v\. To find the energy spectrum and wave func- 
tions, we take 

(33) 
where .E > and U3 and U4 are two-component spinors. 
If M3, U4 are eigenstates of the operator a^h^ — cr^kz 
with eigenvalue k = ^/k"^ + k^ , the energy spectrum is 
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given by E = vik zL g, while if M3, M4 are eigenstates of 
the same operator with eigenvalue —k, the energies are 
given by £^ = —vik ± g. We thus see that the energy 
vanishes, not at the origin {kx,kz) = (0,0), but on the 
circle yfcl-l-fcf = g/vi- This circle of points, along with 
_E = 0, is not invariant under Lorentz transformations 
which will be discussed at the end of this section. Note 
that if V2 had been positive and equal to Vi, the energy 
spectrum of Eq. p2p would have taken a Lorentz invari- 
ant form as discussed below. 

We now consider a wave incident on the line junction 
from region 1. Then there are four possibilities; it can get 
reflected back to region 1 with amplitude r, transmitted 
to region 2 with amplitude i, and transmitted to regions 
3 and 4 with amplitudes t' and t" respectively. Since 
there are four amplitudes to be found and four matching 
conditions in Eq. (IHUl) . the scattering problem can be 
solved. In general, all the four amplitudes will be non- 
zero. However, let us now consider the case of normal 
incidence; then we indeed find that r = f — as expected 
due to the conservation of a^. To be explicit, the incident 
wave here is of the form 



fpi 



1 

71 



Aky-Et) 



(34) 



where E — vik. Conservation of energy equal to E, mo- 
mentum along the x direction equal to zero, and a^ equal 
to 1 in all the regions imply that the transmitted waves 
in regions 3 and 4 must be of the form 



V'3 = ^[t'[ { e-'=^ + t" ; 



V'4 = i[t'(^j)e-'='^-t"(^})e-'="^]e-^*, 

(35) 

where k' = {E/vi) — g and k" — {E/vi) + g, and we have 
normalized the eigenstates corresponding to k' and k" to 
unity. The expressions in Eq. p5|) are obtained by finding 
solutions of Eq. (j32[) with energy E and positive group 
velocities dE/dk' and dE/dk" so that they propagate 
in the —z direction away from the line junction. Using 
Eq. dSOl) along with Eqs. (IM1I55)) gives t' = t" = l/^/2. 
Note that \r\^ + \t\^ -K |i'p + \t"\^+ = 1 as desired. 

The case of a general angle of incidence can be studied 
in a similar way. We assume that the plane waves associ- 
ated with the incident wave in region 1 with amplitude 1 , 
the reflected wave in region 1 with amplitude r, the trans- 
mitted wave in region 2 with amplitude t, and transmit- 
ted waves in regions 3 and 4 with amplitudes t' and t" are 

given bv e'^i^^^+^vV~-E*) ^i(k^x-kyy-Et) ^i(k^x+kyy~Et) 

^i{k^x-k z-Et) g^^^ ^i(k^x-k z-Et) j-gspectively, where 
ky > 0, E 



Vl 



jU2 I 1,2 y 



Wi9-E/v,y 



k" = ^J{g + E/viY — k'^, and we have assumed that 
< E/vi < g and g ± E/vi > \kx\- The above expres- 
sions and conditions ensure that the group velocity of 
the incident wave points toward the line junction but the 



group velocities of the other four waves point away from 
the line junction as desired. We will not write down the 
corresponding two-component wave functions here, but 
simply note that the boundary conditions in Eq. pop 
will provide four equations which will determine the val- 
ues of r, t, t' and t" . In general, all these four amplitudes 
will be non-zero. In the limit that the angle of incidence 
e -^ ±7r/2, we flnd that r ^ -I and t, t', t" -^ 0. 

We may return to the case of surface velocities with the 
same sign as discussed in Sec. II, and ask how the results 
for scattering from the line junction would change in that 
system if one introduces additional surfaces 3 and 4 cou- 
pled by a tunneling amplitude as in Fig. [51 Considering 
Eq. ([5^ with V2 = wi, we find that a tunneling ampli- 
tude g gives rise to an energy gap, instead of producing 
gapless states as in the case of velocities with opposite 
signs. To be specific, the energy spectrum of Eq. ([5^ 
would be given by E"^ = v\{kl. + k1)+g'^, so that \E\ has a 
minimum value of g. (Note that this form remains invari- 
ant under Lorentz transformations based on the velocity 
Vl. For instance, E"^ — v\{k'^ + fc^) remains invariant un- 
der the transformation kx — > {kx — vE/v\)/ ^Jl — v^ lv\, 
E ^ {E — vkx)/y/l — v'^/vf and k^ -^ k^, for any value 
of w whose magnitude is smaller than vi). As a result, an 
electron which is incident from region 1 with an energy 
less than g will either reflect back to region 1 or transmit 
to region 2, and the surfaces 3 and 4 will not play a role 
in this scattering problem. 



V. DISCUSSION 

To summarize, we have shown that a number of inter- 
esting phenomena can occur at a line junction separating 
two surfaces of topological insulators. Assuming invari- 
ance under time reversal, a line junction can have an arbi- 
trary parameter which makes the Hamiltonian Hermitian 
and is consistent with current conservation; this parame- 
ter may be understood as arising from a potential barrier 
lying along the junction. The conductance through the 
line junction depends on this parameter, and there may 
also be edge states which propagate along the junction. 
For a junction which separates two surfaces which are 
at an angle with respect to each other, we have shown 
that the conductance depends on the bending angle, and 
there are also edge states which lie along the junction. 
Our analysis has been based on simple models in which 
the junction width d has been taken to be zero; as a re- 
sult, many properties of the edge states, such as their 
velocity and spin orientation, depend only on dimension- 
less parameters such as a and but not on the magnitude 
of the momentum k along the edge. Realistic junctions 
would have a finite width; the effects of the width need 
to be studied. 

The edge states which lie along the junction differ in 
two ways from the plane wave states which lie far from 
the junction. The velocity of the edge states is always 
less than the velocities of the plane waves, so that the 



energy of the edge states is smaller in magnitude from 
the energy of the plane waves for a given value of the 
momentum along the line junction. Secondly, the orien- 
tation of the spin of the edge states is generally different 
from the orientation of the plane waves; in particular, 
the spins of the edge states have a non-zero component 
in the direction perpendicular to the surfaces. These dif- 
ferences in the energy-momentum dispersion and the spin 
structure imply that it may be possible to observe various 
features of the edge states and to distinguish them from 
the plane waves (which reside far away from the line junc- 
tion) by using spin-resolved angle-resolved photoemission 
spectroscopy and tunneling from a spin polarized STM 
tip placed close to the junction. 

Finally, for the case in which the velocities on the two 
surfaces have opposite signs, we have explicitly demon- 
strated that a wave incident on the junction from either 
side must generally transmit into some additional sur- 
faces which must emerge from the junction. 

Before ending, we would like to point out that our anal- 
ysis has ignored the effects of the hexagonal warping of 
the dispersion which is known to occur in TIs'^^ . It would 
be interesting to study the effects of warping on the con- 
ductance through the junction and on the edge states. 
We would also like to mention recent studies of scatter- 
ing and bound states near ferromagnetic domain walls 
on the surface of a TI^, transport across step junctions 
on a TI surfac o^^'^^ , edge states induced by a magnetic 
field at junctions of two TI surfaces^'?, bulk models for 
studying states at surfaces with arbitrary orientations^^, 
and tachyon-like states at the interface of two TIs— . 
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Appendix A: General boundary condition which 
ensures Hermiticity 

In this Appendix, we will discuss the general linear 
condition between (xi)i/-!.o- and {x2)y^o+ (and similarly 
for (xi)j;-i.o- and {x2)y^o+) which satisfies Eq. ([3]). Let 
us assume that 



iX2)y^0+ = U (Xl)y^o- 

(X2)y^o+ = U (x'i)y^o- 



(Al) 



where U is an arbitrary 2x2 matrix. Then Eq. ^ will 
be satisfied if 



Using the fact that the identity matrix / and the three 
Pauli matrices a form a basis for 2 x 2 matrices, we find 
that the general solution to Eq. (|A2p is given by 



U 



exp[— iacr^ — i/3 — ja"^ — 6a^ 



(A3) 



where a, /3, 7 and S are real parameters; Eq. (jA2[) follows 
from the fact that [/V^ = cr^C/-i. (Eq. ^M^ may be 
contrasted with the form of a unitary matrix V which 
satisfies V^'^V = I; such a matrix can be written in terms 
of four real parameters as V = exp[— iacr^ — i(3 — ija^ — 

We will now present a model of a junction which pro- 
vides a physical understanding of the four parameters 
appearing in Eq. (|A3[) . Consider a thin barrier of width 
d extending from y = to y — d where there is a constant 
term V consisting of both potential and magnetic parts, 



V{y) = A + Ba'' + Ca^ +D(yy, 



(A4) 



where A, B, C and D are real so that the Hamiltonian is 
Hcrmitian. We can think of ^ as a potential, and B, C 
and D as being proportional to the three components of 
a magnetic field which has a Zeeman coupling to the spin 
of the electron. We will eventually be interested in the 
limit of a J-function barrier, so that ^4, B, C, Z? — ?► 00 
and d — >■ keeping Ad, Bd, Cd and Dd fixed. 

The Hamiltonian in the region Q < y < d is now given 
by 

Hd = -iV2 (Cr'^dy - (7^9:,) + A + B(T''+C(7'+D(jy. 

(A5) 
Let us look for a state with energy E and momentum k 
in the x direction; the corresponding wave function is of 
the form 

i^ix,y,t) - fiy) e^^'^^-^'l (A6) 

In the region < y < d, f{y) satisfies the equation 

[-iv2 cr'^dy - V2kay + A + Ba"" + Ca'- + L>cr^] / 



Ef. 



(A7) 



c/V"[/ 



(A2) 



Since we are interested in the limit A, B, C, D ^^ 00, 
we can ignore in Eq. (jA7p the terms of order E and k 
which are finite. The solution for / is then given by 

f{y) = eM{y/v2){-^Acy''-^B-Cay + Da'')]f{Q), 

(A8) 
which implies that 

f{d) ^ cxp[{d/v2)i-iAa'' - iB - Ca^ + Dcr')] /(O). 

(A9) 
We now take the limit A, B, C, D ^ 00 and d ^ 
keeping Ad = U2a, Bd = V2/3, Cd = V2"f and Dd = ~V26 
fixed. By superposing states with different values of E 
and k, we then see that any wave function 4'i^i Vi t) must 
satisfy 

ij{x,0+,t) = exp[-ia(T'' -if3 -ja^ - Sa""] ip{x,0-,t). 

(AlO) 
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This is the general boundary condition proposed in 
Eqs. dH]) and ([Ml) . 

Let us now impose invariance under time reversal. We 
then expect that only the potential A can remain non- 
zero in Eq. (jA5l) while the other three terms must be 
zero since they can be interpreted as arising from a mag- 
netic field. To see this formally, we observe that under 
time reversal, we have to complex conjugate the Haniil- 
tonians in Eq. Q and wave functions, and transform 
the time t — ^ — t; in addition, we have to do a unitary 
transformation by a^ in order to maintain the form of 
the Hamiltonian. In other words, idtip = Hip goes to 
idtaVip* = ayH*ayil)*, where ayH*ay = H. On apply- 
ing this transformation to Eq. (jA5|) . we see that A can 
be non-zero but B, C and D must vanish if V has to be 
time reversal invariant. Hence the parameters /3, 7 and 
5 must vanish in Eq. (JAIOI) . We are therefore left with 
the single parameter a given in Eq. Q. 

If a 7^ 0, the wave functions satisfying Eq. ^ are dis- 
continuous at 2/ = even if V2 — vi. This is not surpris- 
ing. We recall that for the Schrodinger equation which is 
second order in spatial derivatives, a 5-function potential 
barrier leads to a discontinuity in the first derivative of 
the wave function. For the Dirac equation which is first 
order in spatial derivative, a (5-function potential leads to 



a discontinuity in the wave function. 

We would like to emphasize here that our aim has been 
to obtain the most general boundary condition consistent 
with certain Hamiltonians separated by a line junction. 
By doing so, we have included the effects of a thin bar- 
rier which could possibly be present at the junction. In 
a particular system of interest, the relevant parameters 
defining the boundary conditions would have to be de- 
termined from either the microscopic parameters of the 
system or from experimental studies of either transport 
across the junction or the edge states present there. For 
a given system (for instance, if there is actually no junc- 
tion present so that the junction width d is exactly zero), 
it may turn out that all the four boundary parameters 
are equal to zero; however, this would only be a special 
case of our analysis. 

To conclude, our analysis shows that in general four 
real parameters are required to completely specify the 
boundary condition; this generalizes the discussion in 
Refs. l28l - [3ll where no such parameters were considered. 
A similar situation is known to arise for the Schrodinger 
equation with a discontinuity at one point; the general 
boundary condition at that point which ensures Her- 
miticity and conserves the current also has four parame- 
ters as discussed in Refs. yl andliol. 
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